Tests of a High Resolution Time of Flight System 
Based on Long and Narrow Scintillator 
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We have tested a prototype time-of- flight system based on bulk scintillator block of dimensions 2.5 x 2.5 x 200 
cm 3 . Using a calibration scheme similar to the one used in actual collider experiments, we have achieved 
a resolution of 71 ps using Amperex XP2020/UR photomultipliers and 81 ps using proximity-focusing fine- 
mesh photomultipliers (Hamamatsu R2021). Results are also obtained for scintillating fiber blocks of the 
same dimensions. Good internal reflectivity of the bulk scintillator block resulted in resolutions superior to 
the fibre blocks. A single- photon pulsed laser system was used to study photomultipliers and the results were 
used in a Monte Carlo simulation of the system to study the critical elements that determine the resolution. 



1 Introduction 

In many high energy physics experiments the identi- 
fication of particle species is accomplished through 
time-of-fiight (TOF) measurements utilizing a sys- 
tem of plastic scintillators read out by photomulti- 
plier tubes (PMT's). A typical configuration for a 
detector at a colliding beam facilities is a barrel ge- 
ometry consisting of a single layer of long rectangular 
scintillation counters parallel to the beam, segmented 
in azimuth, and read out at both ends by photomul- 
tipliers coupled to light guides [[!]. Two important 
factors that define the performance of such TOF sys- 
tems are the timing resolution and the efficiency to 
register a clean hit. Using 2.8 meter long counters 
of thickness 5 cm and width 10 cm, the CLEO II 
experiment obtains resolutions of 139 ps for bhabha 
events and 154 ps for hadronic events || El. The 
efficiency for obtaining good timing information for 
a track is limited primarily by events where two or 
more particles strike the same counter. In these cases 
one is no longer able to effectively utilize the timing 
measurements from both ends, and the resolution is 
consequently substantially worse. The efficiency at 
CLEO II for spherical events is approximately 85% 
and for jetlike events, approximately 75%. 

The factors that determine the timing resolution 
are quite well understood: the number of photoelec- 
trons detected, ^| the dispersion of path length from 
the scintillation point to the photocathode, the scin- 
tillation decay time, and the transit time spread of 



the photomultiplier. Thus, one tries to use fast and 
efficient scintillator with long attenuation length, fast 
photubes with small transit time spread covering as 
much area as possible at the end of the scintillator 
preferably without any light guides. As a rule of 
thumb, one loses the precious scintillation photons 
by a fraction proportional to the ratio of the cross 
section of the scintillator to the photocathode area. 
The light guide can add dispersion of photon path 
length and further photon loss ]5j . The lack of light 
guides usually requires photomultipliers that can op- 
erate within a high magnetic field (<~1 T). The thick- 
ness of the scintillator should be as thick as possible 
as long as other detector requirements are met |J . 

The efficiency argues for a fine segmentation. The 
effective attenuation of a finely segmented system, 
however, tends to be small due to the increased num- 
ber of internal reflections. It has been suggested that 
resolutions for long counters can be improved by us- 
ing scintillating fibers ||. The advantage of a fiber 
counter is that the resolution at large distances is not 
severely limited by photon path length dispersion and 
that the loss due to internal reflections are small Q| . 
On the other hand, fiber counters have a small trap- 
ping angle which limits the absolute amount of light 
collected by the PMT. The authors of reference || 
compared bulk and fiber counters and concluded that 
the system using the fiber counter has a better res- 
olution at long distances. However, if the effective 
attenuation length of the bulk counter is made com- 
parable to that of the fiber counter, then the bulk 
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counter would give the superior resolution even at 
long lengths. 

In this article, we report a study of a prototype 
time of flight system using long and narrow scintilla- 
tor counters. Such counters could be arranged in a 
two layer system where the inefficiency due to over- 
lapping tracks would be substantially reduced. More- 
over, the timing resolution of such double-layer sys- 
tem would be roughly factor of \/2 better than that of 
a single-layer system (9| . In such system, one also need 
to deal with the technique to install PMT's whose 
diameters are larger than the transverse size of the 
scintillator block, and to combine multiple measure- 
ments for a given track. These are sujects of future 
studies. 



2 Experimental setup 

The three different scintillator counters studied have 
the same dimensions of 2 meters x 2.5 cm x 2.5 cm. 
The first was bulk Bicron BC-408 scintillator (decay 
time 2.1 ns) wrapped in aluminized mylar. The sec- 
ond was a 6 x 6 bundle of 4 mm diameter Bicron 
BCF-10 scintillating fibers (decay time 2.4 ns) with 
a polystyrene-based core and polymethylmethylacry- 
late cladding, glued together with epoxy. The indices 
of refraction were 1.60 for the core and 1.49 for the 
clad. The third counter was similar to the second, ex- 
cept that it was composed of 2 mm diameter fibers. 
For both types of fiber, the clad thickness was 3% 
of the core diameter. The bulk scintillator was ini- 
tially cast between glass planes. The unfinished sur- 
faces were then milled by diamond fly cutting |ll|| . 
This finishing technique yields a reflectivity superior 
to that obtained by standard polishing, achieving an 
improved effective attenuation length of the counter. 
A study of reflectivity of NE110 scintillator |[o| shows 
that any polishing dramatically deteriorates the re- 
flectivity. The ideal fabrication method would be to 
cast all four lengthwise surfaces against glass plates. 

Two different types of PMT's were used to detect 
the scintillation light from the test counter. The first 
was a modification of the popular Amperex XP2020 
called the XP2020/UR. It has a 12 stage dynode and 
a bialkali photocathode of useful diameter 4.5 cm. 
The XP2020/UR differs from the XP2020 in the in- 
ternal connections of accelerating and focusing grids 
and it allows a higher voltage to be applied between 
the cathode and the first dynode to reduce the tran- 
sit time. The base design is shown in Figure [j] |lS^] . 
A further improvement in resolution was obtained by 
adjusting the two potentiometers at Gl and D2 to 
maximize the potential drop from K to Gl and from 
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Figure 1: Base design for the Amperex XP2020/UR. 
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Figure 2: Base design for the fine-mesh photo multi- 
plier Hamamatsu R2021. 



Dl to D2, respectively. As will be shown later in 
the single photoelectron study, the Gl voltage is crit- 
ical in making the transit time uniform across the 
photocathode surface. We applied —3000V to the 
XP2020/UR when it was used with a fiber counter, 
and —2500V with the bulk counter to avoid satura- 
tion. Such a tube would have to be placed outside the 
magnetic field of a detector, and would be coupled to 
the counters via light guides. 

The second type was a Hamamatsu R2021 
proximity-focusing fine-mesh PMT which can be 
placed inside the detector's magnetic field. This PMT 
has a 12 stage dynode structure and a 2.2 cm diame- 
ter bialkali photocathode. The gain is specified to be 
5 x 10 4 at gauss, and 1.5 x 10 3 in a 10 kilogauss field 
parallel to the axis of the PMT. We amplified the out- 
puts of the R2021's with Mini-Circuits ZFL1000-LN 
low noise xlO amplifiers. The base design for this 
PMT is shown in Figure |. When a fine-mesh PMT 
is operated inside a high magnetic field, the primary 
effects of the field are (a) to increase the probability 
of electrons to hit the wire of the mesh due to the 
spiraling motion along the field, (b) reduce the prob- 
ability of secondary emission electrons to reach the 
next stage, and that (c) the cloud of amplified elec- 
trons follows roughly the direction of magnetic field 



2 



as it passes through the layers of meshes JT3]. The 
effect (b) is the main reason for gain loss in a high 
field, while the effect (a) result in higher gains when 
the PMT axis is at some angle 9 (about 30 degrees) 
with respect to the field than zero angle. The ef- 
fect (c) is responsible for the sharp drop of gain for 8 
greater than some value (about 50 degrees) . We have 
not studied fine-mesh PMT's inside high magnetic 
field. Studies show, however, that the deterioration 
in timing resolution is not large |l3|, 

The scintillator and PMT to be tested were placed 
in a cosmic ray telescope as shown in Figure |3|, op- 
tically coupled using an UV-transparent grease GE 
Viscasil 60000. The cosmic ray events were selected 
by coincidence of the trigger counters with dimen- 
sions 2 x 2 x 1.2 cm 3 each. The reference counters 
positioned next to the trigger counters provided a 
measurement of the time that the cosmic ray passed 
through the test counter. All four reference and trig- 
ger counters were mounted on a trolley which could 
be positioned anywhere along the test counter. 

The coincidence of trigger PMT pulses defined the 
TDC start time and generated a gate for the ADC. 
The ADC pedestal was measured every five minutes. 
Each pulse from the two reference PMT's and the 
two test PMT's was split and passed to both a LeCroy 
2249A ADC and leading edge discriminators (LeCroy 
623B and 4608C) whose output was sent to a LeCroy 
2228A TDC calibrated to 50 ps/count. The discrimi- 
nator threshold for the reference counters was chosen 
to be —30 mV, which is roughly 5% of the typical 
pulse height, in order to maximize acceptance. The 
threshold for the trigger counters was chosen to be 
—400 mV, which was roughly 1/3 of the typical pulse 
height in order to exclude background events. 

3 Attenuation length 

Figure ^ shows the pulse height distributions of the 
test counter PMT's. The trolley was positioned at 
the center. Results are shown for four different com- 
binations of PMT's and counters: the XP2020/UR's 
at -2500V on the bulk counter, the R2021s on the 
bulk counter, the XP2020/UR's at -3000V on the 4 
mm diameter fiber counter, and the XP2020/UR's at 
—3000V on the 2 mm fiber counter. 

The effective attenuation length of each counter 
was measured by fitting a exponential function to 
the pulse height as a function of the distance from 
the each end of the counter. This gives two indepen- 
dent measurements for a single scan. Figure |^ shows 
the results obtained for the bulk counter, indicating 
an effective attenuation length of 160 cm. For the 
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Figure 4: Distribution of ADC counts for a test 
counter photo multiplier (PMT3) for four differ- 
ent configurations: the XP2020/UR's on the bulk 
counter, the R2021s on the bulk counter, the 
XP2020/UR's on the 4mm fiber counter, and the 
XP2020/UR's on the 2mm fiber counter. 

fiber counters, we find an effective attenuation length 
168 cm for the 4 mm diameter fiber counter, and 147 
cm for the 2 mm diameter fiber counter. Note that 
the effective attenuation length of the bulk counter 
is comparable to the effective attenuation lengths of 
the fiber counters. Thus, relative to the fiber coun- 
ters, the bulk counter does not suffer from the loss of 
photon statistics at the locations far from PMT. As 
a result, the bulk counter may have an overall timing 
resolution superior to those of the fiber counters due 
to the better photon statistics, which we observe to 
be the case as reported below. 



4 Timing resolution 

In the following discussion, we assign the numbers 1 
and 2 to the reference counter PMT's, 3 to the test 
counter PMT at z = cm, and 4 to the test counter 
PMT at z = 100 cm. The mean time measured by 
the two reference counters is the time at which the 
cosmic ray crossed the test counter. The time mea- 
sured by PMT i must be corrected to account for the 
time required for the pulse to reach the discriminator 
threshold (the 'time- walk' correction). The correc- 
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Figure 3: Cosmic ray telescope test setup. 

We collect data with the trolley placed at 25 cm in- 
tervals along the test counter, from z = 25 cm to z 
= 175 cm. We combine the data from all z positions 
and determine a single value for a\ and for a 2 which 
is independent of z. 

Using the values for a\ and 02, we solve for the 
constants a 3 , a 4 , c 3 and C4 of the test counter PMT's 
by minimizing the widths of the distributions 

t l (z) = T3 - c 3 z - r rof+ 

and 

T R (z) = T4 — c 4 z - T ref+ 

with 

T re{+ = (Ti + T^)/2. 

The constant a is the inverse of the effective speed 
of light in the test counter and z is the distance to 
the PMT in question. Note in particular that the 
constants 03 and 04 are independent of z. We thus 
follow the method of calibration identical to the one 
used at the CLEO II detector where the beam cross- 
ing time is replaced by T ro f + . We do not require that 
C3 = c 4 following the standard calibration technique 
of CLEO II. The value of the time- walk constant were 
chosen to be 63 = 64 = 0.15 as described later. We 
define the time resolutions of the left and right PMT's 
to be 

= V^M-^cf 

where <r^(z) and <Jb,{ z ) are the measured widths of 
Tl(z) and Tr(z), respectively. We then construct 
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Figure 5: Attenuation length of the bulk counter, 
determined by fitting the mean pulse height from each 
end to an exponential. 

tion has the form 

T! = Ti- ai Q7 b * 

where Ti is the measured TDC time and Qi is the 
pedestal-subtracted ADC count. We use 61 = 6 2 = 
0.5 for the reference PMT's. The constants a\ and 
a<i are determined from the data by minimizing the 
width of the distribution 

T ref _ = (T[ - 



4 



I 200 



I 000 



800 - 



600 - 



400 - 



200 



Ts ys te m at z = 1 00 cm 



a = 1 .876 ± 0.015 counts 
= 93.8 ± 0.8 ps 



.... 

520 




530 540 
TDC counts 



550 



Figure 6: Distribution of T sys t e m for the trolley at z = 
100 cm, with the XP2020/UR's on the bulk counter. 
The reference counter resolution must be subtracted 
in quadrature from the width of this distribution to 
obtain <r to{ (z = 100). 



the weighted average of crossing time measurements 
as follows: 



T. 



system 
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The measured width of this distribution is denoted 
as tT sys tom(z)- Figure || shows T systC m for the bulk 
counter at z = 100 cm, using XP2020/UR's. The 
distribution is well-described by a gaussian, and has 
a resolution of 94 ± 1 ps. 

Finally the contribution from the reference coun- 
ters must be subtracted in quadrature to obtain the 
resolution of the test counter: 



Otof(z) = vA 



system 



In previous studies of similar time-of-flight systems, 
the resolution of the reference counter time difference 
T re f_ was often taken to equal the average time T re f + . 
However, such an assumption can be made only if 
the cosmic ray tracks have a narrow angular devia- 
tion and in the limit of small trigger counters. With 
our geometry, the angular acceptance of cosmic rays 
is approximately 20 degrees from vertical. This gives 
rise to a correlation between the two reference timings 
that cancels in T lc f + but not in T re f_ since when the 
track is close to one of the PMT's then it is away from 
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Figure 7: Reference counter configurations. The 
standard configuration (left), and that for the mea- 
surement of reference counter resolution (right). The 
effect of timing correlation between Reference 1 and 
2 due to the track path for Tl + T2 in the left con- 
figuration is the same as that for Tl - T2 in the right. 



the other (see Figure 0). Under these conditions, the 
resolution of T re f- is significantly greater than the 
resolution of T re f+. There is also a smaller correla- 
tion introduced by the variation of the track position 
in the mid plane between the two trigger counters. In 
order to properly account for these correlations, the 
reference resolution was measured by placing PMT's 
facing opposite directions. The width of T re f_ in this 
configuration should be equal to that of T re f + when 
both PMT's face the same direction. This assump- 
tion was verified by Monte Carlo studies fl6fl . The 
reference time resolution was thus measured to be 
<7 re f = 59 ± 1 ps (Figure ||). Using this value of a lc f , 
we obtain the resolution of the test counter system 
to be CTtof = 73 ps at z = 100 cm. The contribution 
to atof from the finite TDC bin size is negligible. 

Having defined the timing resolution of the sys- 
tem, we examined the resolution obtained for differ- 
ent thresholds on the test counter discriminators. For 
this study, data were collected at z = 100 cm. The 
test counter discriminator thresholds were set to a 
fixed percentage of the average pulseheight for each 
end, and runs were taken for various thresholds in the 
range of 1% to 40% of the average pulse height. Fig- 
ure |9j shows the resolution versus discriminator level 
for the bulk and the 4 mm fiber counters. The thresh- 
olds were chosen to be 5% for the 4mm fiber counter, 
and 10% for the bulk counter. 

We next determined the value of the time- walk con- 
stants 63 and 64 constraining them to be equal. For 
this study data were collected for all z positions, and 
the simple average of Ctof (z) for all z position was cal- 
culated. Figure |lC| shows how the average resolution 
depends on 63 for the bulk and for the 4 mm fiber 
counter, each at their optimal discriminator thresh- 
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Figure 8: Resolution of reference counters, defined as 
the width of T ref - = (T{ - T%)/2 when the reference 
counter PMT's are oriented in opposite directions. 
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Figure 10: 0t o f(z = 100) as a function of the expo- 
nential term of the test counter time-walk correction 
63 (= 64): bulk counter (right axis) and 4mm diame- 
ter fiber counter (left axis). 
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Figure 9: 0t o f(z = 100) as a function of discrimina- 
tor threshold for test counter PMT's: bulk counter 
(solid square) and 4mm diameter fiber counter (open 
triangle). 



olds. The data favor lower values than the oft-used 
value of 0.5. We choose 63 = b 4 = 0.15 for both fiber 
and bulk counters using the XP2020/UR's. 

Figure |ll| shows 03(2), 174(2), and atoi(z) for the 
bulk counter using the XP2020/UR's. The left 
counter resolution was found to be worse than that 
of the right counter. This was found to be consis- 
tent with the photon statistics as measured and at- 
tributed to an imperfect optical coupling of the PMT 
to the scintillator block. We take the simple average 
of <r to f (z) to obtain an overall resolution of 71 ± 1 ps 
for the bulk counter. Figures 12 and [l3] show simi- 
lar data for the XP2020/UR's on the 4 mm diameter 
and 2 mm diameter fiber counters, respectively. We 
find an average otaf °f 125 ± 3 ps for the 4 mm fiber 
counter and 133 ± 2 ps for the 2 mm fiber counters. 

Lastly, we studied the bulk counter using Hama- 
matsu R2021 proximity- focus fine- mesh PMT's. For 
this study, we used the discriminator threshold at 
10% of the pulse height. The optimal value for 
the constants 63 and 64 were found to be 0.30 us- 
ing the method described earlier. Figure |lj shows 
03(2), (J 4( Z ): an d <Jtof(z) for the R2021s on the bulk 
counter. The average a to f is found to be 81 ± 1 ps 
which is slightly worse than the resolution obtained 
by XP2020/UR. It should be noted, however, that 
the photocathode of R2021 covered only 65% of the 
end of the bulk scintillator and the difference is con- 
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Figure 11: Time resolutions of the left PMT (0-3(2)), 
the right PMT (04(2:)) and the weighted average 
(crtof(z)). For the bulk counter with XP2020/UR's. 
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Figure 13: Time resolutions of the left PMT (0-3(2)), 
the right PMT (0-4(2)) and the weighted average 
(otof(z)). For the 2 mm diameter fiber counter with 
XP2020/UR's. 
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Figure 12: Time resolutions of the left PMT (o 3 (z)), 
the right PMT (0-4(2)) and the weighted average 
(otof(z)). For the 4 mm diameter fiber counter with 
XP2020/UR's. 
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Figure 14: Time resolutions of the left PMT (o 3 (z)), 
the right PMT (0-4(2;)) and the weighted average 
(atof(z)). For the bulk counter with R2021s. 
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Transit Time vs. XP2020/UR Photocathode Position 



Configuration 


Resolution (ps) 


XP2020/UR on bulk 


71 ±1 


XP2020/UR on 4 mm fiber 


125 ±3 


XP2020/UR on 2 mm fiber 


133 ±2 


R2021 on bulk 


81 ±1 



Tabic 1: The measured time resolution <7t f for each 
configuration averaged over the length of the test 
counter. 
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Figure 16: Transit time dependence on the location 
of the incident photon on the photocathode position 
for XP2020/UR. 



Figure 15: The setup for the single photon study. 
The entire system except the TDC logic is placed in 
a dark room. 

sistent with the photons lost. Table |] summarizes the 
average Otof for each configuration. 

5 Single photon study 

Examination of the single photon response of the 
XP2020/UR is instrumental in understanding one of 
the components of the time-of-flight system resolu- 
tion. In particular, it is useful to know how the tim- 
ing resolution varies according to the position on the 
photocathode at which the photon is incident. In ad- 
dition, it also allows the gain of the PMT/base com- 
bination to be determined accurately. These are used 
as input parameters for the Monte Carlo simulation 
of the system. 

The light source used for this investigation was a 
Laser Photonics LN120C nitrogen laser, which deliv- 
ered pulses of wavelength 337.1 nm and r.m.s. dura- 
tion 70 ps, each with 70 /iJ, and at a repetition rate 
of about 10 Hz. After exiting the laser, the beam was 
split (see Figure |l5|). Figure |l5|. One beam triggered 
a fast photodiode (Newport 818 BB, with rise time < 
200 ps) and defined the start time of the event while 
the other beam proceeded through a series of neutral 



density filters and then to the PMT. The attenua- 
tion provided by the filters was sufficiently large that 
a PMT pulse resulted only once for every ~10 laser 
pulses. Thus, approximately 95% of the signal events 
were due to one photoelectron. 

The PMT was mounted on a movable x-y table 
which allowed the position of the incident photon to 
be varied across the photocathode. The PMT was 
magnetically shielded by mu metal inside degaussed 
soft iron. An iris with a 1 mm diameter aperture was 
placed in front of the PMT to define the beam size 
at the photocathode. The signal from the PMT was 
sent to a Mini-Circuits ZFL-1000LN amplifier and 
then split to an ADC and TDC. In order to shield 
the PMT from electrical noise generated by the laser, 
each of the two apparati was separately enclosed in an 
electrically grounded box of 1/8 inch-thick aluminum. 

The PMT was oriented such that longer dimension 
of the first dynode was horizontal. The repositioning 
of the stage was done in darkness to keep the PMT 
on and stable for the duration of the scan. A common 
pulse height correction was applied to all runs from 
the scans of both axes. 

Figure |l^ shows the mean TDC stop time for hor- 
izontal and vertical scans of the photocathode. We 
have observed that the shape of this curve depends 
on the potential applied to the accelerating grid Gl. 
Figure [l7| shows the mean corrected TDC stop time 
as a function of horizontal beam location for three 
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Figure 17: Transit time vs the location of incident 
photon on photocathode for XP2020/UR at -3000V. 
Plotted for three different potentials of the accelerat- 
ing grid Gl. 



different Gl potentials. The position dependence of 
the transit time is clearly minimized for large poten- 
tial differences between photocathode and Gl. The 
variation of the shape with Gl potential suggests that 
the potential can be adjusted to provide a nearly uni- 
form transit time for all photocathode positions, and 
that the resolution can be degraded if this parameter 
is not optimized. 

The gain of each PMT/base combination was de- 
termined from the mean number of ADC counts mea- 
sured per pulse, accounting for the 0.25 pC/count cal- 
ibration of the ADC and for the splitter. The gain 
was measured at two different HV values, and found 
to vary as (HV) n with n between 8.5 and 9.7. This is 
consistent with the dependence expected for a twelve 
stage PMT. 



6 Monte carlo simulation 

We developed a Monte Carlo program to simulate 
the performance of the system with the XP2020/UR 
PMT's on the bulk counter in order to understand 
the relative importance of the physical processes con- 
tributing to the resolution of the system. 

Cosmic ray tracks were generated with a cos 2 {aff) 
angular distribution where 6 is the polar angle with 
respect to the vertical. The coefficient a = 2 was cho- 
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Figure 18: Monte Carlo effective attenuation length 
of bulk counter as a function of the reflectivity. The 
best estimate of reflectivity was taken to be 0.9925. 



sen to match the shapes of the reference counter pulse 
height as described later. The number of photons 
was generated by the Landau distribution. The pho- 
ton emission sites were distributed randomly along 
the length of the track, and the emission directions 
were isotropic. The decay time of the scintillator was 
taken to be 2.1 ns. 

The path of each photon was traced through the 
scintillator until the photon either reached the PMT, 
or until it was lost by one of a number of different 
mechanisms. Photons were absorbed in the scintilla- 
tor assuming a bulk attenuation length of 3.8 meters 
for BC-408. The index of refraction of BC-408 was 
taken to be 1.57 and the reflection/transmission at 
boundaries were modeled according to Fresnel's law 
for unpolarized light. Furthermore, the photon was 
assumed to be lost at the boundary with probabil- 
ity 1 — R due to imperfection. The photomultiplicr 
window was modeled as a circular piece of glass with 
index of refraction 1.47 in direct contact with the 
scintillator. The quantum efficiency of the photo- 
cathode was assumed to be 26% averaged over the 
photon emission spectrum. 

The reflectivity R was chosen as follows: The effec- 
tive attenuation length of the bulk counter is a func- 
tion of the reflectivity, the bulk attenuation length, 
and the geometry of the counter. Figure [l8] shows the 
effective attenuation length as a function of the reflec- 
tivity We chose R = 0.9925, which is the value that 
results in the measured effective attenuation length 
of 1.58 meters. 
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Figure 19: Comparison of Monte Carlo and data for 
the pulse height distributions of the reference counter 
(Ql) and the test counter (Q3). The horizontal scales 
are adjusted by roughly 20% to compare the shapes. 

Two transit time effects were modeled. First, the 
transit time difference between center and edge of 
the photocathode was assumed to vary as r 2 . For 
the reference counter, we took this difference to be 
250 ps. For the XP2020/UR, we took the results 
from Figure [l6|. Second, the transit time spread was 
modeled by a gaussian of sigma 250 ps. 

The distribution of photoelectron arrival times was 
convoluted with a function which described the time 
response of the PMT to a delta function light in- 
put. For the test counter PMT's, this function was 
the pulse shape produced by the XP2020/UR on the 
single photon apparatus, measured by a 4 GHz Tek- 
tronix SCD5000 transient digitizer. For the reference 
counters, this function was a gaussian of width 1.5 
ns. 

The absolute yield of photons per unit path length 
was determined by requiring that the Monte Carlo 
pulse height distribution matches with the data dis- 
tributions. We find a production rate of 5890 pho- 
tons/cm. This corresponds to a BC-408 light output 
which is 20% of anthracene. Using the resulting pulse 
shape for each track crossing, the discriminator trig- 
gering time was defined to be the time the pulse shape 
crosses the threshold. 

The pulse height distributions of the reference 
counter (Ql) and the test counter (Q3) are shown for 
Monte Carlo and data in Figure [l{j. The lower shoul- 
der for the reference counter corresponds to tracks 
with large dip angle 9 which intersect the sides of the 
reference counters. The coefficient a in the angular 
distribution cos 2 (a6) was determined to match the 
shape and size of this shoulder. 

We calibrate the simulated data in the same way as 
described previously. We find a to i — 55 ± 3 ps for the 



(a) # 7 /cm % of 


Otof 


anthracene (ps) 


1440 5 


93 ±6 


2880 10 


72 ±6 


5880 (nominal) 20 


55 ±3 


8650 29 


37 ±7 


11540 38 


32 ±9 


23080 77 


33 ± 5 


(b) TTS (ps) 


0tof (ps) 


1000 


91 ±3 


500 


61 ±4 


250 (nominal) 


55 ±3 


125 


53 ±5 





52 ±5 


(c) Scintillation 


Otof 


decay time (ps) 


(PS) 


2.5 


55 ±6 


2.1 (nominal) 


55 ±3 


1.75 


51 ±6 


1.5 


48 ±6 


1.0 


46 ±7 


0.0 


40 ±7 


(d) Center/Edge 


Otof 


difference (ps) 


(PS) 


400 


56 ±4 


300 


55 ±4 


200 (nominal) 


55 ±3 


100 


57 ±4 





55 ±4 



Table 2: Monte Carlo study: Variation of the param- 
eters in the simulation and the resulting change in 
otof = 100). The parameters varied are: (a) num- 
ber of photons produced per unit path length of track 
also expressed as a percentage of anthracene output, 
(b) the XP2020/UR transit time spread, (c) the BC- 
408 decay time, (d) the XP2020/UR photocathode 
center/edge transit time difference. 
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z = 100 cm trolley position, in reasonable agreement 
with the data. 

We varied a number of parameters in the Monte 
Carlo and observed the resulting changes in atoi(z — 
100). The parameters included the number of pho- 
tons produced per unit length of cosmic ray track, 
the scintillator decay time, the transit time spread, 
and the photocathode center-edge difference of the 
transit time. In all cases, the variations were ap- 
plied to the test counter and its PMT's, but not to 
the reference counters. Tables ^ summarizes these 
results. We find that the resolution is limited pri- 
marily by photon statistics. In our application, the 
resolution is fairly insensitive to the transit time dif- 
ference between photocathode center and edge. This 
is because only the central portion of the photocath- 
ode of XP2020/UR is used. In applications that use 
the entire photocathode, it will be critical to adjust 
the potential of Gl to make the transit time uniform 
across the photocathode surface. 



7 Conclusions 

We have measured the resolution of a small cross sec- 
tion prototype time-of-flight system. We find that a 
bulk counter manufactured to give an effective atten- 
uation length comparable to a fiber counter gives a 
better resolution than the fiber counter due to the 
improved photon statistics. Results using magnetic- 
field-resistant photomultiplier tubes are comparable 
to those obtained with a fast photomultiplier of a 
more conventional design. We have studied the sin- 
gle photon response of XP2020/UR photomultiplier 
and found that the grid voltage adjustment is critical 
in flattening the distribution of transit time as a func- 
tion of photon position on photocathode. This is par- 
ticularly important if the entire photocathode is to 
be illuminated. We developed a detailed monte carlo 
simulation program which reasonably reproduces the 
data. The monte carlo shows that the critical element 
is the number of photoelectrons, emphasizing in par- 
ticular the importance of the internal reflectivity of 
the scintillator block. 
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